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Lithography provides the precision to pattern large arrays of metallic nanostructures with varying 
geometries, enabling systematic studies and discoveries of new phenomena in plasmonics. However, surface 
plasmon resonances experience more damping in lithographically-defined structures than in chemically- 
synthesized nanoparticles of comparable geometries. Grain boundaries, surface roughness, substrate effects, 
and adhesion layers have been reported as causes of plasmon damping, but it is difficult to isolate these 
effects. Using monochromated electron energy-loss spectroscopy (EELS) and numerical analysis, we 
demonstrate an experimental technique that allows the study of these effects individually, to significantly 
reduce the plasmon damping in lithographically-defined structures. We introduce a method of 
encapsulated annealing that preserves the shape of polycrystalline gold nanostructures, while their 
grain-boundary density is reduced. We demonstrate enhanced Q-factors in lithographically-defined 
nanostructures, with intrinsic damping that matches the theoretical Drude damping limit. 



W 



hy do plasmonic structures have poorer performance when they are fabricated with top-down tech- 
niques compared to those fabricated with bottom-up approaches? 



Surface plasmon resonances produce strong electric fields within sub-diffraction-limited volumes that can be 
applied in various nanophotonic devices. Using top-down lithographic approaches such as electron-beam 
lithography (EBL), these resonators can be patterned deterministically over large arrays with sub- 10 nanometer 
dimensional control 1 " 3 . As a universal damping parameter for resonators, the quality-factor (Q-factor) has been 
extensively studied with different plasmonic nanostructures 4,5 . Unfortunately, in addition to damping due to the 
intrinsic properties of the material, extrinsic sources such as grain boundaries, surface roughness, and adhesion 
layers 6 " 10 also introduce damping and suppress Q-factors in lithographically-defined nanostructures. In contrast, 
wet-chemically- synthesized metal nanoparticles can be monocrystalline and atomically smooth 1112 and possess 
higher Q-factors. However, due to the challenges in precise placement of these particles, lithography provides an 
acceptable compromise for plasmonic device applications. To enhance Q-factors of lithographically-defined 
structures, thermal annealing 71314 , and substituting metal adhesion layers with organic ligands 915 have shown 
promise. However, more detailed studies are needed, e.g. to anneal-out grain boundaries without shape change, 
and suppress Rayleigh instabilities that break up long and narrow nanostructures 1617 . As high Q-factor plasmonic 
resonators are useful in practical applications such as surface- enhanced Raman scattering (SERS) 18 , color print- 
ing 19 , and light extraction from quantum emitters 20 , efforts to improve Q-factors are relevant. 

In this work, we present a readily- applicable method 21 that combines the design freedom of top-down 
lithography and the high Q-factors of chemically-synthesized nanoparticles, where the intrinsic damping is at 
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the Drude limit 711 ' 22 . By incorporating a thermally-stable but easily 
removable encapsulation layer, the EBL-defined nanostructures can 
be annealed to reduce grain boundaries without changing the geo- 
metry or surface roughness of the nanostructures. We applied this 
method to isolate and analyze the contributions of the various damp- 
ing sources in plasmonic nanostructures. Ensemble- averaging 
effects were avoided in our analysis by studying single nanostructures 
with monochromated electron energy-loss spectroscopy (EELS) in a 
transmission electron microscopy (TEM) 23 " 28 . Interestingly, our 
results show that despite the incomplete removal of grain bound- 
aries, annealed nanostructures have damping parameters close to 
those of monocrystalline nanoparticles, matching the Drude damp- 
ing parameter of bulk metal. We also show the importance of a 
hitherto unexamined parameter, i.e. the contact area between the 
nanostructure and a dielectric substrate. 

Results 

Figure 1 shows the process flow for reducing the grain boundary 
density of nanostructures. Au nanostructures with a thickness of 
25 nm were fabricated directly onto 30 nm thick, free-standing 
Si 3 N 4 TEM support membranes 29 ' 30 . A 30 nm thick layer of hydrogen 
silsesquioxane (HSQ, XR-1541 2% solids) was then spin-coated 
onto the samples. Thermal annealing was performed in a Jipelek 
rapid thermal processor under nitrogen environment at 400 °C for 
10 minutes. After annealing, the HSQ encapsulation was removed by 
immersing the sample in a 2% vol. solution of hydrofluoric acid (HF) 
for 15 s. The shape and surface roughness of the nanostructures were 
preserved due to the encapsulation. In contrast, control samples 
without encapsulation experienced changes in shape, size, and sur- 
face roughness 16 as can be seen in Fig. SI of the Supplementary 
Information (SI). 

The nanostructures were fabricated onto Si 3 N 4 membranes to 
enable the direct high- resolution imaging of the grain boundaries. 
Figure 2 shows TEM images of ring- and rod-shaped nanostructures 
with a 1 nm Cr adhesion layer, before and after annealing. The 
highly-polycrystalline nature of the gold nanostructures prior to 
annealing is seen in the TEM image of Fig. 2(a) with grain sizes of 
— 10 nm. Grain boundaries are marked by dashed, yellow lines in the 
close-up panels. As shown in Fig. 2(b), the same region after anneal- 
ing had fewer grain boundaries, likely due to grain-boundary migra- 
tion and Ostwald ripening at elevated temperatures 31 . The average 
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Figure 1 | Process flow for reducing grain boundaries of 
lithographically-defined nanostructures with minimal change in their 
shape and surface roughness, (a) Gold nanostructures on 30 nm thick, 
amorphous Si 3 N 4 membranes, (b) Encapsulation by a 30 nm thick layer of 
HSQ. (c) Annealing at 400°C for 10 minutes, (d) Immersing the sample in 
dilute HF to remove the encapsulating HSQ. 



diameter of the grains increased to —20 nm, leading to a significant 
reduction of the total number of crystals, as evidenced by the reduced 
number of bright spots in the electron diffraction patterns (bottom 
panels, Fig. 2(a) and (b)). 

Shape preservation was afforded by HSQ encapsulation and not by 
adhesion layer pinning. An example of a structure fabricated without 
Cr is shown in Fig. 2(c). The length of this nanorod was 454 nm 
before and 446 nm after annealing, a reduction in length of —2%. 
The edge profiles of the nanorod were still preserved by the encap- 
sulation, confirming the minimal effect of the process on the shape 
and surface roughness of the nanostructures. 

EELS measurements were done on the same structures both before 
and after they were annealed; the experimental procedures are 
described in the SI. EELS results for the nanorod of Fig. 2(c) are 
presented in Fig. 3(a). Up to four peaks corresponding to different 
resonant harmonic modes are evident, and their corresponding plas- 
mon maps are shown beneath the respective peaks. The effects of 
annealing are evident in the enhanced intensity and narrowed full- 
width at half maximum (FWHM) of the resonant peaks. Q-factors 
were extracted from the spectra by using a previously- described 
method 32 . The blue shift of the plasmon peaks observed in Fig. 3(a) 
after annealing is likely caused by a combination of factors, including 
the reduction of the rod length during annealing, and possibly the 
thinning of the substrate during etching, as is shown in more detail in 
Fig. S3 in the SI. 

Numerical calculations using FDTD were performed to extract the 
Q-factors and field distributions of the nanostructures, as described 
in the Methods section. The relative permittivity for bulk Au was 
described using the Drude-Lorentz model 7 , expressed as follows: 



or 



1 + i(xT p 



(i) 



Here co p is the plasma frequency, T p is the damping parameter for 
the Drude term, co m , and T m are the resonant frequency and the 
damping parameter of the m th Lorentz oscillator, respectively. The 
increase in damping due to factors such as grain boundaries and 
impurities is accounted for by an increase in the Drude damping 
parameter. Following the work of Liu et at. 10 , a damping factor a 
was introduced to quantify this increase relative to bulk Au, where 
a = 1 is defined to be at the Drude limit. By fitting experimental 
permittivity data 33 with two Lorentzian terms, we obtained the 
following parameters: co p = 9 eV, a = 1, Y p = 0.07 eV, fi = 
0.3, Wi = 2.7 eV, I~\ = 0.3 eV, f 2 = 0.8, co 2 = 3.05 eV and T 2 
= 0.5 eV, as shown in Fig. S2 in the SI. Figure 3(b) shows spectra 
for the fundamental mode of a 100 nm long gold rod with a = 1, 
1.3, and 1.8. Q-factors were extracted from fitted Lorentzian func- 
tions to the simulated spectra 4 . 

Experimental Q-factors extracted from three sets of gold nano- 
structures are shown in Fig. 3(c) as a function of resonance energy. 
To verify the effects of the adhesion layer and grain boundaries on the 
damping process, we compared these experimental results with 
simulated Q-factors of nanostructures with lengths from 40 to 
275 nm and square cross sections 25 nm wide. Structures with Cr 
adhesion have the lowest Q-factors, and could be fitted with a damp- 
ing factor of a = 1.8. Without Cr, Q-factors immediately improved 
by —25%. Notably, samples with Cr adhesion did not show a signifi- 
cant increase in Q-factor (<10%) after annealing, though it was 
evident that they experienced a reduction in grain boundary densi- 
ties, as seen in Fig. 2. It is likely that Cr diffuses into the gold during 
annealing, to form Cr-rich crystallites 34 , creating new scattering sites 
that compensate the reduction in damping from grain boundaries. 
Conversely, the samples without Cr exhibited an improvement in Q- 
factor by as much as —40% after annealing (Fig. 3(c) inset). The as- 
fabricated structures without Cr could be fitted with a = 1.3, which 
reduced to a = 1 after annealing. 
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Figure 2 | Reducing the density of grain boundaries with shape preservation of EBL-defined nanostructures. Panels (a, b) show structures with Cr 
adhesion layers, while the rods in panel (c) have no adhesion layer. TEM image and electron diffraction pattern of the same ring structure before (a) and 
after annealing (b), showing a reduction in the density of grains and grain boundaries, (c) TEM images of a rod structure before and after annealing. 



Discussion 

Annealing strongly enhances the Q-factors for plasmons with res- 
onance energies below —1.7 eV. This effect can be understood from 
the relatively low- frequency Drude damping parameter in Au, Y p ~ 
0.1 eV. Therefore, only resonances that occur at low frequencies are 
most strongly affected by a modification to this parameter. 
Consistently, Fig. S2 in the SI illustrates that the damping factor 
affects only the imaginary part of the permittivity at energies below 
— 1.7 eV. The decrease in Q-factors > 1 .7 eV is due to damping from 
interband transitions in Au 4 . 

Figure 4(a) shows that the wet- chemically synthesized nanorods 
have comparable but systematically higher Q-factors than the 
annealed, EBL-patterned nanostructures. We identified two main 
reasons for this observation: (1) smaller physical cross sections and 
volumes of the rounded single-crystal nanorods 11 ' 35,36 , and (2) the 
smaller contact area of these nanorods with the substrate. The widths 
of the measured wet-chemically synthesized nanorods were in the 
range of 17 to 32 nm while the EBL patterned nanostructures were 
systematically wider, ranging from 23 to 32 nm. To enhance mutual 
comparison, the simulations used 25 nm as the lateral widths for all 
the simulated structures. As shown in Fig. 4(a), numerical simula- 
tions using a = 1 for structures with a pentagonal cross section (blue 
line) agree well with experimental data for the single-crystal nanor- 
ods. In reality, the cross sections of our wet- chemically synthesized 
particles are neither pentagonal nor perfectly circular 37 , but some- 
where in-between, see for example Fig. S4 in the SI. As a result, their 
Q-factors would fall between the blue and black solid lines in 
Fig. 4(b), the latter of which simulates nanostructures with rounded 
corners. 

The normalized E-field magnitude distributions of gold nano- 
structures with two different lengths for points I-IV are simulated 
in the panels on the right of Fig. 4. For resonances around 1 eV 
(panels I, II), the field distributions are comparable for both the 
square and pentagonal cross sections. At these low energies, radiative 



damping dominates. As a consequence, their relatively smaller 
volume leads to a somewhat higher Q-factor for pentagonal nanor- 
ods 32 . This effect is confirmed in Fig. 4(b), where at low energies, 
structures with smaller cross-sectional areas (i.e. smaller volumes) 
have higher Q-factors. Conversely, for high-energy resonances 
around 2.1 eV (panels III, IV), absorptive damping dominates and 
the stronger field localization at the metal-substrate interface dic- 
tates that structures with the smallest contact areas have the highest 
Q-factors, as seen in Fig. 4(b). 

The effect of the substrate is further studied by comparing Q- 
factors for structures with and without substrates, also plotted in 
Fig. 4(b). For a given cross section, the presence of a substrate 
improves the Q-factor at low energies, but strongly reduces the Q- 
factor at high energies. The different effects of the substrate arise 
again from the interplay between radiative and absorptive damping. 
The presence of a substrate reduces radiative damping at low ener- 
gies, by localizing fields at the metal-substrate interface, and leads to 
an overall increase in Q-factor, see also Figures S5 and S6 in the SI. 
On the other hand, this substrate effect leads to a decrease in Q-factor 
at high energies due to the enhanced absorptive damping at the 
metal-substrate interface. The sensitive dependence on the metal- 
substrate contact area can also explain the scatter in the data of 
Fig. 4(a) for single- crystal nanorods: Some of these rods are likely 
to contact the substrate only at localized points, unlike the conformal 
contact in lithographically-defined structures. 

It is worth pointing out that even though grain boundaries remain 
after annealing, the damping factors of the annealed structures still 
match those of single-crystal nanorods. We predict that the effect 
still of grain boundaries on damping strongly depends on their ori- 
entation, their position along the particle length with respect to the 
plasmon resonance nodes, and whether low- or high-angle grain 
boundaries are present 38 . The crystallographic defect density can 
vary greatly between different types of grain boundaries. For an 
isotropic material such as gold, defect-free twin boundaries 
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Figure 3 | Effect of encapsulated annealing and the adhesion layer on plasmon damping, (a) Experimental EELS spectra before and after annealing for 
the single nanostructure shown in Fig. 2c. The spectra were first normalized at the 'zero-loss peak' maximum, to allow mutual comparison of their 
intensities. Peaks I-IV with their corresponding EELS plasmon maps underneath show the first four harmonic resonance modes, (b) Simulated extinction 
efficiency showing peak broadening with increasing damping factor a. (c) Q-factors measured from over 20 nanostructures from samples fabricated with 
a Cr adhesion layer (green triangles), without any adhesion layer before annealing (blue squares), and after encapsulated annealing (red disks). Solid lines 
are numerical simulations with a chosen to achieve good agreement with experimental data. Inset: Q-factor enhancement as a function of resonance 
energy. 



potentially introduce less damping than high-angle grain boundaries 
(having grain misorientations >10°) loaded with crystallographic 
defects that act as scattering centers for electrons 31 . We therefore 
hypothesize that the large plasmon damping in as-deposited samples 
is predominantly due to the presence of high-angle grain boundaries. 
As grain-boundary mobility increases with the degree of misorienta- 
tion, these high-angle grain boundaries are also most easily removed 
by annealing. Our results suggest that the boundaries that remain are 
predominantly low-angled ones that have little effect on plasmon 
damping. The annealed structures, while polycrystalline, are there- 
fore optically indistinguishable from a single crystal. Finally, over the 
entire energy range that was investigated, surface roughness of the 
nanostructures was found to have very little effect on the Q-factor 4,7 , 
as can be seen in Fig. S7 in the SI. 

In answer to our introductory question, we demonstrated that 
lithographically- defined nanostructures experience more plasmon 
damping than wet- chemically synthesized metal nanoparticles due 
to the combined effect of grain boundaries, the larger metal-sub- 
strate interface, and metal adhesion layers. We presented a facile 
process to improve the Q-factors of nanostructures that are fabri- 
cated using top-down lithography. Encapsulation by HSQ, followed 
by annealing and HSQ removal resulted in polycrystalline gold nano- 
structures with Q-factors comparable to those of single- crystal 
nanorods. Both structures could be modelled using damping para- 
meters of bulk gold. Our results provide design guidelines that indi- 
cate the dominant source(s) of damping for a given energy range, and 
will enable the fabrication of high Q-factor plasmonic resonators 
and antennas. If left intact, the HSQ encapsulation could also prevent 



degradation of metal resonators operating at elevated temperatures, 
such as in near-field transducers for heat-assisted magnetic record- 
ing, and photoelectron-emission cathodes. 

Methods 

EELS experiments were performed in scanning TEM (STEM) mode using an FEI 
Titan TEM with a Schottky electron source. The microscope was operated at 80 kV, 
using a STEM probe with a diameter around 1 nm. With a Wien-type monochro- 
mator, an energy resolution with typical FWHM values of 70 meV was obtained. 
Mapping and spectroscopy was done with a Gatan Tridiem ER EELS detector, 
applying a 7 mrad collection semiangle using STEM-based spectrum imaging and a 
modified binned gain averaging acquisition routine 39 . The EELS spectra in Fig. 3 were 
taken by summing all EELS spectra from the spectrum image in the first 5 nm off the 
surface of the whole nanostructure. The resulting summed EELS spectra are therefore 
representative of the nanostructure as a whole. The spectra were normalized at the 
'zero-loss peak' maximum, making it possible to mutually compare the spectral 
intensity before and after annealing. The background signal was measured from a 
bare amorphous Si 3 N 4 TEM support membrane, fitted to- and subtracted from the 
experimental EELS plasmon spectra. EELS maps show the loss signal integrated over 
an energy window of 0.1 eV energy around the plasmon peaks of interest. Q-factors 
were calculated from the first and— when possible— higher- order plasmon modes 40 in 
the EELS spectra, as described in detail elsewhere 32 , by applying Olivero and 
Longbothum's accurate empirical approach 41 . Q-factors were extracted from the 
most prominent peaks for a given nanostructure, up to the third-order harmonic. 

Synthesis of gold nanorods. Gold(III) chloride trihydrate (HAuCl 4 -3H 2 0, >99.9% 
trace metals basis), Sodium citrate dihydrate (>99%), hexadecyltrimethylammonium 
bromide (CTAB, >99%, powder), and L-ascorbic acid (Biotra, >99.0%) were 
purchased from Sigma- Aldrich Singapore Pte Ltd. Ultra-pure deionized water 
(ELGA pure lab water systems) was used to prepare all aqueous solutions. The 
preparation of seeds —3.5 nm in diameter was done following the three- step seed- 
mediated method introduced by Jana et al 42 . A 20 ml of gold seeds aqueous solution 
containing HAuCl 4 -3H 2 0 (2.5 X 10" 4 M) and tri-sodium citrate (2.5 X 10" 4 M) was 
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Figure 4 | Effect of the substrate and the interface on plasmon damping, (a) Experimental Q-factor plot for annealed, EBL-defined structures (red 
disks), and chemically synthesized, single- crystal nanorods (blue pentagons). The right-hand panels show the I El field distributions, normalized at their 
maxima (side view and edge-on views (inset)) for points I-IV. (b) Simulated Q-factors for structures with a square cross section (red lines), a pentagonal 
cross section (blue line), and a rounded cross section (black line), all using a = 1. Dashed lines are simulated Q-factors in the absence of the Si 3 N 4 
substrate. 



prepared in a 22 ml glass vial. Then, 0.6 ml of an aqueous ice-cold NaBH 4 (0.1 M) 
was added at once into the growth solution. The solution was allowed to stir for 
2 hours before proceeding to the next step. A 30 ml of growth solution containing 
CTAB (0.1 M) and HAuCl 4 *3H 2 0 (2.5 X 10" 4 M) was prepared. This growth 
solution was distributed into three polypropylene tubes (labeled A, B and C) with 
9 ml each. Next, 0.05 ml of ascorbic acid (0.1 M) was added into each tube followed 
by gentle inversion mixing. Then, 1.0 ml of the gold seeds solution was added into 
tube A. After 15 s, 1.0 ml from tube A was transferred into tube B. After another 15 s, 
1.0 ml from tube B was transferred from tube B into tube C. Solution in tube C which 
contained gold nanorods, spheres and plates, was subjected to centrifugation for 
shape separation. Typically, 1.5 ml of the particle solution was centrifuged at 
2000 rpm for 6 min repeatedly. The supernatant was poured off and the solid part 
containing nanorods and plates was re-dispersed in about 0.1 ml of deionized water, 
from which a droplet was dried on a Si 3 N 4 TEM support membrane. 

Simulations. The commercial finite-difference time-domain (FDTD) code from 
LUMERICAL was used to simulate the damping parameters and field distributions. A 
plane wave with the electric field polarized along the long axis of the nanorod 
illuminates a nanorod at normal incident angle. Perfectly matched layer (PML) 
boundary conditions encompass one isolated nanorod. For samples fabricated by 
EBL, a square cross section, 25 nm wide, was set for all the relevant simulations 
regardless of annealing condition or adhesion layer, while for the wet-chemically 
synthesized particles, a 25 nm wide pentagonal cross section was applied. Electric 
field plots were normalized to their maxima. Q-factors were obtained by fitting the 
fundamental mode of the extinction spectra with a Lorentzian function, using 
nanorods with lengths varying from 40 to 275 nm. The relative permittivity for Si 3 N 4 
was set to be 4. 
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